Localization from quantum interference in one-dimensional disordered potentials 
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We show that the tails of the asymptotic density distribution of a quantum wave packet that 
localizes in the the presence of random or quasiperiodic disorder can be described by the diagonal 
term of the projection over the eigenstates of the disordered potential. This is equivalent of as- 
suming a phase randomization of the off-diagonal/interference terms. We demonstrate these results 
through numerical calculations of the dynamics of ultracold atoms in the one-dimensional speckle 
and quasiperiodic potentials used in the recent experiments that lead to the observation of Anderson 
localization for matter waves [Billy et at, Nature 453, 891 (2008); Roati et at, Nature 453, 895 
(2008)]. For the quasiperiodic case, we also discuss the implications of using continues or discrete 
models. 
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Since the seminal work by P. W. Anderson the 
behavior of quantum particles in disordered potentials 
has been a central topic in physics @, Q. In particu- 
lar, the diffusion of a quantum wave packet in a one- 
dimensional disordered systems has been longly inves- 
tigated from the theoretical point of view ^-12], and 
also observed experimentally in two recent experiments 
with matter waves [l^ . In these experiments, the 
initial wave packet is prepared as the ground state of a 
trapped Bose-Einstein condensate. Then, the trapping is 
switched off and the condensate is let expand in a one- 
dimensional wave guide in the presence of disorder, in a 
regime where interatomic interactions can be neglected. 
Under certain conditions, disorder induced exponential 
localization (Anderson localization) is observed. Similar 
results have been reported for t he p ropagation of light in 
quasiperiodic photonic lattices [15| . 

Usually, the localization is interpreted in terms of lo- 
calization of single plane wave components, as a result 
of multiple scattering from the random potential ^ J^] . 
Recently, a more elaborate approach has been discussed 
in jl^ . In this article we propose an alternative and 
complementary viewpoint, in terms of the eigenstates of 
the disordered potential V{x). We show that the tails 
of the disorder averaged asymptotic density distribution 
are determined by the diagonal term of the projection of 
the initial wave packet on the eigenstates of the single 
particle Hamiltonian H — /2m + V{x), and that this 
is equivalent to an effective randomization of the eigen- 
state phases. In fact, though the interference term does 
not vanish (actually, its amplitude is always of the or- 
der of the diagonal term), it contributes just with local 
fluctuations. 

We demonstrate these results through numerical work 
for the one-dimensional speckle and quasiperiodic po- 
tentials employed in recent experiments with ultracold 
atoms [13;, 14] . This analysis highlights the role of quan- 
tum interference in the localization process, demonstrat- 
ing that it plays a key role for both correlated disorder 



and quasiperiodic potentials, even in the deep localizing 
regime. For the quasiperiodic case, we also discuss the 
different features of continuos and discrete models. 

Model. Let us consider a one-dimensional non- 
interacting system, namely a particle of mass m in 
a potential V{x), described by the hamiltonian H = 
— (?i^/2to)V^ -I- V{x). For simplicity we consider a finite 
size system of length L, whose discrete spectrum is char- 
acterized by eigenvalues En and eigenvectors (l)n{x) — 
{x\n). Then, the evolution of an initial wave packet 
{x\ip{0)) can be expressed in terms of the eigenstates |n) 
as 



(xm)) 
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with c„ = {n\ip{0)). By writing the projection coefficients 
in the modulus-phase representation, c„ = |c„| exp(i0„), 
we can rewrite Eq. ^ as {x\4'(t)) = J2n \cn\4>nix)e^'^"^*\ 
with ipn(t) = On — Ent/Ti. These equations are equivalent 
to solving the Schrodinger equation ihdtipix) = Hipix), 
but have the advantage that the density profile at any 
time can be computed directly on the knowledge of just 
the initial conditions. The density distribution n[x,t) = 
|(a;|V'(t))P at any time t is given by 

n 

+ ^ |c„||c„.|0:(a;)0„'(x)e^^'^"'.'(*) 

= ridiag (x) + Hoff {x, t) (2) 

where we have defined the phase difference A(^„„/ (t) = 
(pn{t) — ^n' [t), and explicitly separated the diagonal and 
the off-diagonal interference terms. The evolution of the 
latter is determined by the phase difference A(^„„'(t), 
whereas the diagonal term does not depend on time. 

Conceptually, in the case of disorder, this approach is 
complementary to the usual description of localization 
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in terms of the localization of single momentum compo- 
nents (plane waves), as a result of multiple scattering 
from the random potential [Hi- Its advantages are that 
each component evolves independently from the others 
(instead, plane waves are not independent - even in the 
non-interacting case - because they are coupled by the 
potential, see also {TS]), and that the time evolution of 
each component can be simply obtained from the knowl- 
edge of the energy spectrum En . From the practical point 
of view, this allows to compute the wave-function at ar- 
bitrary long times with the effort of a single time-step. 

In the limit L ^ oo the spectrum may become con- 
tinuous, or remain discrete, depending on the potential 
V{x). In the first case, it is known that the initial wave 
packet will spread indefinitely, whereas this cannot hap- 
pen if the spectrum is discrete [11]. The discreteness of 
the spectrum can derive both from the presence of a con- 
fining potential, and in this case the boundedness of the 
dynamics is a trivial consequence, or more interestingly 
due to the presence of disorder [2|. The implications of 
this latter case will be discussed in the following. 

Asymptotic regime. Let us consider the effect of the 
interference term in Eq. ([5]), in the asymptotic (long 
times) regime. Initially, the phase relations between dif- 
ferent components are exactly those that reconstruct the 
assigned initial profile. In case of a real wave packet, 
0„ = or TT. Then, they evolve in a deterministic way, 
according to ipnit) = — Ent/h, determining the wave 
packet profile as a result of the interference of many com- 
ponents. For a disordered system in the localized regime 
- characterized by a point (discrete) spectrum - we expect 
the asymptotic distribution, noo{x) = limt^oo n{x,t), to 
be stationary (except for local fluctuations) and deter- 
mined by the localization properties of the eigenstates 
i- 

We argue here that, for a disordered system in the lo- 
calized regime, the behavior of the tails of the asymptotic 
density distribution is determined just by diagonal terms 
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and can be therefore computed by the knowledge of only 
the initial coefficients c„ and the localization properties 
of the eigenstates 4>n{x). In fact, it turns out that the in- 
terference (off-diagonal) term asymptotically contributes 
just with fluctuations around the profile of Eq. ([3]). This 
can be explained in terms of an effective randomization of 
the phases ipn, as will be discussed later on. The approx- 
imate equivalence in Eq. ([3]) refers to local fluctuations 
related to the specific disorder realization, the relation 
becoming more and more accurate when averaged over 
many disorder realizations. 

To characterize the relaxation of the the full density 
on the diagonal term during the evolution of the system, 
we introduce the following quantity 



that measures the average relative weight of the diago- 
nal term (in logarithmic scale). This definition is espe- 
cially useful when the density is exponentially localized. 
Asymptotically, we expect R{t) to approach 1 if ridiagix) 
is a good approximation oi n{x,t). 

As far as concerns the shape of the asymptotic density 
distribution, when the eigenstates are localized around 
Xn — {4>n\x\(f>n) and decay exponentially away from the 
localization center, |0„(x)p ~ exp[— (a: — x„)|/^„], the 
long distance behavior of the asymptotic wave packet 
is determined by the largest localization length i — 
maxlin}. However, the overall profile may not be char- 
acterized by a simple exponential decay, due to the su- 
perposition of many energy components with different 
localization lengths (H, [l2| . 

In the following we illustrate these ideas for the lo- 
calization of matter waves, as considered in recent ex- 
periments [U, [13] • In both cases the spectrum is com- 
puted by mapping the stationary Schrodinger equation 
H(f)n — En4>n ou a discrctizcd grid [l^. All the results 
presented in this paper are averaged over 100 realizations 
of disorder. 

Quasiperiodic potential. Let us start by discussing 
the case of a quasiperiodic potential as considered in 
[S [m (3 • The quasiperiodic potential is created by su- 
perimposing two optical lattices, of wavevector ki and in- 
tensity Vi {i = 1, 2), V{x) = Vi sin^(A;ia;) + V2 sin^(fc2a; + 
9), with Vi 3> V2, corresponding to a primary optical 
lattice weakly modulated by the second one, 6 being an 
arbitrary phase. In this case, it is convenient to ex- 
press lengths in units of the lattice spacing d — ir/ki 
and energies in terms of the associated recoil energy 
Eri = h^kf/2m. In the tight binding regime, Vi S> Em, 
the system Hamiltonian can be mapped onto the discrete 
Aubry-Andre model [3 for the particle amplitude (pj 
at site j of the primary lattice H = —J{(pj_^.iip* + 

c.c.) -I- A cos(27r^j -|- 9)\(pj\'^, where (3 = k2/ki is the ratio 
between the two lattice wavevectors, and the tunneling 
constant J and disorder amplitude A are functions of the 
amphtudes l^i [Ull^]. When the two wavevectors are in- 
commensurate, the model is characterized by a transition 
from extended to localized states, at A/J = 2 [l^. Here 
we chose 13 = (\/5- l)/2, the inverse of the golden mean 
21|, and A/J = 7, in order to discuss a specific example 
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where localization occurs. 

As initial state we consider a gaussian wave-packet 
modulated by the primary lattice, prepared as the non- 
interacting ground state of the primary lattice (V2 = 0) 
plus an additional harmonic trapping (here we consider 
ujko — 0.03 i?_Ri/?i), as in the experiments [l3]. This 
choice guarantees that the initial state \tp{0)) has mainly 
components in first Bloch band of the primary lattice. In 
addition, in order to make a direct comparison with the 
solution of the discrete AA model [^ - that is intrinsically 
a single band model - we also consider a truncated initial 
state ]V'(0)), obtained by removing the components lying 
in higher Bloch bands. 

Let us start from the case of the truncated initial state 
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FIG. 1. (Color online) (a),(b),(c) Asymptotic density distri- 
bution in the quasiperiodic potential, in the localized regime. 
The full expression ^ (red, solid line in (a)) at £ = 10^ Ti/Eri 
is compared with the expression in Q (orange, (b)), and with 
the same Eq. but in case of a random choice of tiie pliases 
(fin (blue line, (c)). All the data refer to an initial wavepacket 
in the single-hand approximation, also shown in (a) (black dot- 
ted line). The straight (dashed) lines in all these three panels 
are obtained from an exponential fit of the density tails in 
(a), with a function oc exp[— 2|a;|/^;oc] (here lioc ~ O.Sld, see 
text), (d) Evolution of the quantity R{t) in Eq. Q as a func- 
tion of time; the dotted horizontal line represents the case of 
Eq. ((2} for a random choice of the phases tpn {R — 0.993). 



IV'(O)). In Fig. [TJa) we show the initial wave packet 
no{x) = |(a;|?^(0))p and the average asymptotic density 
distribution obtained from Eq. ^ at t — 10^ h/ Ej^j. 
The tails of the asymptotic density are characterized 
by an exponential decay, as found from the solution 
of the discrete AA model In particular, by fitting 
the density with a function oc exp[— 2|a;|/£ioc], we find 
iioc ~ 0.81 d, that is consistent with the expected value 
(aa = d/log(A/2J)|A/j=7 ~ 0.8 d for the AA model 
[isj . Then, in Fig. [Tib) we plot the diagonal term 
ndiagix), and in Fig. [TJc) again the full expression of 
Eq. but in case of a random choice of the phases ipn- 
In both cases we keep the same fitting lines of Fig. [Tfa) 
for comparison. These figures show that both approxi- 
mations in (b) and (c) nicely reproduce the actual den- 
sity distribution, confirming that Eq. ([3]), or equivalently 
the phase randomization ansatz, correctly describes the 
asymptotic behavior of the tails in the localized regime. 
Indeed, the three densities in (a), (b) and (c) would 
be hardly distinguishable if drawn together on the same 
plot. From Figs. H^a) and (b) it is easy to evince that the 
ofF-diagonal interference terms are essential in determin- 
ing the evolution of the initial wave packet towards the 
asymptotic density distribution, but once the system is 
in the asymptotic regime they can be disregarded, as long 
as the behaviour of the tails is concerned [l^. The ap- 
proach to the asymptotic regime can be tracked by means 
of the quantity R{t) in Eq. @ , whose evolution is shown 
in Fig. [IJd). It is characterized by a smooth crossover 



FIG. 2. (Color online) Same as Fig. [T] but for the full ini- 
tial wave packet (a;|^(0)), with components also in excited 
Bloch bands (black dotted line, in (a)), (a) Asymptotic den- 
sity distribution from the full expression ^ (red, solid line); 
(b) expression in ([Sjl (orange); (c) Eq. ^ but in case of a 
random choice of the phases ifin (blue). The exponential fit 
~ exp[— 2|a;|/£AA] (dashed line) is the one obtained from Fig. 
ma), (d) Evolution of the quantity R{t) in Eq. Q as a func- 
tion of time, compared with the case of Eq. ^ for a random 
choice of the phases ip„ (dotted line, R ~ 0.971). 



towards i? ~ 1 around a "critical" time ta = 10^ h/Em, 
signaling that for t > ta the s yste m has almost relaxed 
to the asymptotic distribution |23{ . Indeed, the behavior 
shown in Figs. [Ila)-(c) for t = 10^ %/ Em is found at any 
time in the asymptotic regime. Remarkably, we find the 
same behavior even for a single realization. 

In order to discuss the effect of higher bands, that are 
naturally included in the continuous description of real 
experiments, we now consider the full initial wave packet 
{x\il){Q)) [i^- The corresponding asymptotic distribution 
are shown in Figs. [2ja)-(c) (for comparison we also show 
the same fitting lines of Fig. [TJa); note the different scales 
between the two figures). With respect to the previous 
case of a truncated initial state |'0(O)) , there is now a con- 
tinuous background due to the fact that the eigenstates 
in higher Bloch bands are extended. Nevertheless, even 
in this case the full solution is correctly described by the 
asymptotic expression We note that the presence 
of this background represents an intrinsic cut-off on the 
extension of the exponential tails, that is inherent in the 
continuous description of the system. Fig. [2jd) shows 
that in this case the behaviour of R{t), though similar 
to the previous one, presents some quantitative differ- 
ences: initially R{t) starts from a lower value, and then 
saturates slightly below the single band case. Both these 
effects can be attributed to the presence of the extended 
states. In particular, the relatively low initial value is due 
to the very different shape of the initial and asymptotic 
density distributions. 

Speckle potential. Let us now turn to the case of an 
optical speckle potential, that is a correlated disorder po- 
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FIG. 3. (Color online) Localization in the speckle potential, 
(a) Initial wave packet (black dotted line) and symptotic den- 
sity distribution obtained from the full expression in Eq. ((2]) 
(red, solid line), with the tails fitted as oc |a:;|"e~^'''^' (black, 
dashed line); the same lines are plotted also in panels (b) and 
(c) for comparison, (b) Contribution of the diagonal terms, 
see Eq. ([Sjl (orange), (c) Eq. ((2)1 with a random choice of the 
phases </p„ (blue, dashed line), (d) Evolution of the quantity 
K(i) in Eq. Q, compared with the case of Eq. ([2]) for a 
random choice of the phases <^„ (dotted line, R ~ 0.959). 



tential characterized by a negative-exponential distribu- 
tion of the intensity -P(Vb) = e~^"/<^o V(^o), and a spa- 
tial (auto) correlation function r(a:') = 1 -|-sinc(0.88x/^)^, 
^ being the correlation length [2a, [23| • It can be written 
as V{x) = V[)v{x) with the distribution of intensities of 
v{x) being normalized to unit standard deviation. The 
speckles are characterized by the presence of an effective 
mobility edge 'S] , below which the eigenstates localize ex- 
ponentially. Lengths and energies can be conveniently ex- 
pressed in units of the natural scales provided by the cor- 
relation length f and its associated energy — ?i^/2m^^ 
(28j . The speckles can be generated numerically by fol- 
lowing the method described in [13, . 

Here we consider the quantum evolution of an initial 
Thomas-Fermi (TF) wave packet with the macroscopi- 
cally occupied components lying below the effective mo- 
bility edge, as considered in 0, 0, [H, [l^. We choose 
Vb = 0.19 and an initial TF radius i?TF = 340 ^ [H. 
The corresponding asymptotic density distribution and 
the evolution of the quantity R{t) are shown in Fig. |31 
Again, in (a) we plot the asymptotic density as obtained 



from the full expression in Eq. ([2]), that has to be com- 
pared with the expression in Eq. ([3]) containing just the 
diagonal terms, shown in (b), and with the expression 
in Eq. ^ but with a random choice of the phases, in 
(c) . We find that the tails of the density distribution can 
be nicely fitted with the function tiqIx) cx |a;|"e~^'''l^l, as 
discussed in (here we find a 2, 7 « 4 • 10 ^). The 
same fitting lines are reported in all the three panels (a) , 
(b) and (c), for comparison. These figures confirm that 
the phase randomization anstaz holds also for the case 
of the speckle potential, and that the asymptotic density 
distribution can be reproduced by the diagonal term in 
Eq. ([3]). Correspondingly, also the evolution of R{t) in 
Fig. [3Ud) shows a behavior similar to the quasiperiodic 
case. Finally, we notice that a similar behavior holds also 
in case of a single speckle realization, but with larger lo- 
cal density fluctuations with respect to the average case. 

Conclusions. We have shown that a convenient way 
to deal with the quantum localization of a wave packet 
in the presence of disorder is to consider its projection 
over the eigenstates of the disordered potential. With 
this approach, the tails of the asymptotic density distri- 
bution are determined by the diagonal term of the den- 
sity, that is independent of time [30|. This is equivalent 
of assuming a randomization of the phase of the eigen- 
states, and corresponds to the fact that the off-diagonal 
(interference) terms contribute just with local fluctua- 
tions around the diagonal term. As speciflc examples, 
we have considered the cases of two recent experiments 
where Anderson localization was observed for ultracold 
atoms in speckle [11] and quasiperiodic potentials [l4| . 
In the latter case we have also shown that the exponen- 
tial decay of the tails in the continuum description is 
intrinsically limited by the presence of extended states in 
excited Bloch bands. The mechanism discussed in this 
paper provides a complementary approach to the usual 
description in terms of the localization of single momen- 
tum components as a result of multiple scattering from 
the random potential. 
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